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Abstract. Motivated by recent experimental results we reconsider semileptonic D — > PfrVi and 
• D — > V£V( decays within a model which combines heavy quark symmetry and properties of the 

chiral Lagrangian. Using limits of soft collinear effective theory and heavy quark effective theory 
we parametrize the semileptonic form factors. We include excited charm meson states in our 
Lagrangians and determine their impact on the charm meson semileptonic form factors. Then we 
calculate branching ratios for all D — > P£v$ and D — > VlV/ decays. 



The knowledge of the form factors which describe the weak heavy — > light semilep- 
tonic transitions is very important for the accurate determination of the CKM parameters 
from the experimentally measured exclusive decay rates. Usually, the attention has been 
devoted to B decays and the determination of the phase of the V u j, CKM matrix element. 
At the same time in the charm sector, the most accurate determination of the size of 
V cs and V cc l matrix elements is not from a direct measurement, mainly due to theoretical 
Q_i! uncertainties in the calculations of the relevant form factors' shapes. 

Recently, there have been new interesting results on D-meson semileptonic decays. 
The CLEO and FOCUS collaborations have studied semileptonic decays D° — > n~£ + v 
and D° — > K~£ + v JEH]. Their data provide new information on the D° — > n~£ + v and 
D° — > K £ + v form factors. Usually in D semileptonic decays a simple pole parametriza- 
^ . tion was used in the past. The results of Refs. tl|,|2|] for the single pole parameters re- 

quired by the fit of their data, however, suggest pole masses, which are inconsistent with 
the physical masses of the lowest lying charm meson resonances. In their anlyses they 
also utilized a modified pole fit as suggested in yD and their results indeed suggest the 
existence of contributions beyond the lowest lying charm meson resonances III]. 

In addition to these results new experimental studies of charm meson resonances have 
provided a lot of new information on the charm sector BJ, |5l |y, tZD which we can now 
apply to D and D s semileptonic decays. 
The purpose of our studies 

aa 

is to accommodate contributions of the newly discov- 
ered and theoretically predicted charm mesons in form factors which are parametrized 
using constraints coming from heavy quark effective theory (HQET) limit for the region 
of cfinax and in the q 2 ~ region using results of soft collinear effective theory (SCET). 
We restrain our discussion to the leading chiral and l/mn terms in the expansion. 

The standard decomposition of the current matrix elements relevant to semileptonic 
decays between a heavy pseudoscalar meson state \H(pjj)} with momentum p v H and a 
light pseudoscalar meson state \P(pp)) with momentum pp is in terms of two scalar 



functions of the exchanged momentum squared q 2 = (pn — pp) 2 - the form factors 
F + (q 2 ) and Fo(q 2 ). Here F + denotes the vector form factor and it is dominated by 
vector meson resonances, while Fq denotes the scalar form factor and is expected to 
be dominated by scalar meson resonance exchange Hill . In order that the matrix 
elements are finite at q 2 = 0, the form factors must also satisfy the relation ^+(0) = 
F (0). 

The transition of \H(pn)) to light vector meson \ V(py, £y)) with momentum p v and 
polarization vector £y is similarly parameterized in terms of four form factors V, Aq, A\ 
and A2, again functions of the exchanged momentum squared q 2 = (pn — py) 2 . Here 
V denotes the vector form factor and is expected to be dominated by vector meson 
resonance exchange, the axial A\ and A2 form factors are expected to be dominated by 
axial resonances, while Ao denotes the pseudoscalar form factor and is expected to be 
dominated by pseudoscalar meson resonance exchange ifTTIl . As in previous case in order 
that the matrix elements are finite at q 2 = 0, the form factors must also satisfy the well 
known relation A (0) +Ai(0)(m H + mv)/2m v —A2(0)(m H — m v )/2m v = 0. 

Next we follow the analysis of Ref. yfl, where the F + form factor in H — > P transitions 
is given as a sum of two pole contributions, while the Fq form factor is written as 
a single pole. This parametrization includes all known properties of form factors at 
large m#. Using a relation which connects the form factors within large energy release 
approach [12] the authors in Ref. J3] propose the following form factor parametrization 

(1 — x) { 1 — ax) l—bx 

where x = q 2 / . 

Utilizing the same approach we propose a general parametrization of the heavy to light 
vector form factors, which also takes into account all the known scaling and resonance 
properties of the form factors. As already mentioned, there exist the well known HQET 
scaling laws in the limit of zero recoil H 1 311 while in the SCET limit q 2 — ► one obtains 
that all four H — > V form factors can be related to only two universal SCET scaling 



functions [12] 



The starting point is the vector form factor V, which is dominated by the pole at 
t = mjj* when considering the part of the phase space that is close to the zero recoil. For 
the heavy — > light transitions this situation is expected to be realized near the zero recoil 
where also the HQET scaling applies. On the other hand, in the region of large recoils, 
SCET dictates the scaling described in [12]. In the full analogy with the discussion 
made in Refs. ylllj], the vector form factor consequently receives contributions from 
two poles and can be written as 

Vtf) = n C " ( '- g) . . (2) 
(1 — x) { 1 — ax) 

where x = q 2 /m 2 ^* ensures, that the form factor is dominated by the physical H* pole, 
while a measures the contribution of higher states which are parametrized by another 
effective pole at m 2 ^ = m 2 H */a. 

An interesting and useful feature one gets from the SCET is the relation between V 
and Ai lfl2LH5lllfl. 12] at q 2 ~ 0. When combined with our result 0, it imposes a single 



pole structure on A\. We can thus continue in the same line of argument and write 



1-b'x 



(3) 



Here % = mjj / '{mn + my) 2 is the proportionality factor between A\ and V from the 
SCET relation, while b' measures the contribution of resonant states with spin-parity 
assignment 1 + which are parametrized by the effective pole at m?,,» = mjj^/b'. It can 

eff 

be readily checked that also A\, when parametrized in this way, satisfies all the scaling 
constraints. 

Next we parametrize the Aq form factor, which is completely independent of all 
the others so far as it is dominated by the pseudoscalar pole and is proportional to a 
different universal function in SCET. To satisfy both HQET and SCET scaling laws we 
parametrize it as 



again parametrizes the contribution of higher pseudoscalar states by an effective pole at 
m H' = m H I a '• T ne resem blance to V is obvious and due to the same kind of analysis |Q] 
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although the parameters appearing in the two form factors are completely unrelated. 

Finally for the A2 form factor, due to the pole behavior of the A\ form factor on one 
hand and different HQET scaling at g^ ax on the other hand, we have to go beyond a 
simple pole formulation. Thus we impose 



which again satisfies all constraints. Due to the relations between the form factors we 
only gain one parameter in this formulation, b" . This however causes the contribution of 
the 1 + resonances to be shared between the two effective poles in this form factor. 

At the end we have parametrized the four H — > V vector form factors in terms of the 
six parameters c' H , a, a', b', c" H and b". 

In our heavy meson chiral theory (HM^T) calculations we use the leading order 
heavy meson chiral Lagrangian in which we include additional charm meson resonances. 
The details of this framework are given in [8] and Ja]. We first calculate values of the 
form factors in the small recoil region. The presence of charm meson resonances in 
our Lagrangian affects the values of the form factors at g^ ax and induces saturation 
of the second poles in the parameterizations of the F + (q 2 ), V(q 2 ) and Ao(q 2 ) form 
factors by the next radial excitations of D/\ and mesons respectively. Although 
the D mesons mat not be considered heavy enough, we employ these parameterizations 
with model matching conditions at gj^ax- Using HQET parameterization of the current 
matrix elements |8l|9[], which is especially suitable for HM^T calculations of the form 
factors near zero recoil, we are able to extract consistently the contributions of individual 
resonances from our Lagrangian to the various D — > P and D — > V form factors. We use 
physical pole masses of excited state charmed mesons in the extrapolation, giving for the 




Ai{q 2 ) 



(ma + my)E,c / H ( \ — a) +2myc'J J ( \ — a') 
(m H -m v )(l-b'x)(l-b"x) 



(5) 



pole parameters a = mjj t /mjj /t , a' = mjj/mfj,, b' = m^/mjj . Although in the general 
parameterization of the form factors the extra poles in F + , V and Aq i 2 parametrized 
all the neglected higher resonances beyond the ground state heavy meson spin doublets 
(0~, 1 _ ), we are here saturating those by a single nearest resonance. The single pole 
q 2 behavior of the A\(q 2 ) form factor is explained by the presence of a single 1 + state 
relevant to each decay, while in A2(q 2 ) in addition to these states one might also account 
for their next radial excitations. However, due to the lack of data on their presence we 
assume their masses being much higher than the first 1 + states and we neglect their 
effects, setting effectively b" = 0. 

The values of the new model parameters appearing in D — > Plvi decay amplitudes ii 
are determined by fitting the model predictions to known experimental values of branch- 
ing ratios -> K~i+v), &(D+ -> K°£+v), ^(D° -> 7T£ + v), &(D+ -> 7t°£+v), 

In our calculations of decay widths we 
neglect the lepton mass, so the form factor Fq, which is proportional to q^, does not 
contribute. For the decay width we then use the integral formula proposed in [19] with 
the flavor mixing parametrization of the weak current defined in Ja]. 

Similarly in the case of D — > VlV\ transitions we have to fix additional model pa- 
rameters [9] and we again use known experimental values of branching ratios SS{Dq — > 
K*-£+v), &(D+ -> 4>£+v), &(D+ -> p°£+v), &(D+ -> K*°£+v), as well as partial 
decay width ratios T L /T T (D + -> K*°£+v) and T + /r_(D+ -> K*°£ + v) [ 1 8 ] . We cal- 
culate the decay rates for polarized final light vector mesons using helicity amplitudes 
H+ - o as in for example [20]. By neglecting the lepton masses we again arrive at the in- 
tegral expressions from 11911 with the flavor mixing parametrization of the weak current 
defined in [9]. 

We first draw the q 2 dependence of the F + and Fq form factors for the D° — > K , 
D° — ► n~ and D s — > K° transitions. The results are depicted in Fig.[U Our model results, 
when extrapolated with the double pole parameterization, agree well with previous 
theoretical ElLIZal and experimental studies whereas the single pole extrapolation 
does not give satisfactory results. Note that without the scalar resonance, one only gets 
a soft pion contribution to the Fq form factor. This gives for the q 2 dependence of Fq a 
constant value for all transitions, which largely disagrees with lattice QCD results 11221 
as well as heavily violates known form factor relations. 

We also calculate the branching ratios for all the relevant D — > P semileptonic decays 
and compare the predictions of our model with experimental data from PDG. The results 
are summarized in Table [TJ For comparison we also include the results for the rates 
obtained with our approach for F + (q 2 nSLX ) but using a single pole fit. It is very interesting 
that our model extrapolated with a double pole gives branching ratios for D — > P£V£ in 
rather good agreement with experimental results for the already measured decay rates. 
It is also obvious that the single pole fit gives the rates up to a factor of two larger than 
the experimental results. Only for decays to r\ and 77 ' as given in Table[H an agreement 
with experiment of the double pole version of the model is not better but worse than for 
the single pole case. 

We next draw the q 2 dependence of all the form factors for the D — > K~*, D — > p~ 
and D s — > <j> transitions. The results are depicted in Fig.EJ Our extrapolated results for 
the shapes of the D — ► V semileptonic form factors agree well with existing theoretical 



studies [20, 2JJ, |22l |24], while currently no experimental determination of the form 



TABLE 1. The branching ratios for the D — ► P semileptonic decays. Comparison of our model fit 
with experiment as explained in the text. 



Decay 


x-r^i /- ill 11 i\r ^t/ n 

S3 (model, double pole) [%] 


/ 11* 1 l\r f~rf ~1 

SS (model, single pole) [%] 


J? (Exp. [18]) [%] 


D° — » K 


3.4 


4.9 


3.43 ±0.14 




0.27 


0.56 


0.36 ±0.06 


D+-TJ 


1.7 


2.5 


2.5±0.7 


D+-Tj' 


0.61 


0.74 


0.89 ±0.33 


D+^K° 


9.4 


12.4 


6.8±0.8 


D+->7Z° 


0.33 


0.70 


0.31 ±0.15 


D+ -> TJ 


0.10 


0.15 


<0.5 




0.016 


0.019 


< 1.1 




0.20 


0.32 





factors' shapes in these decays exists. 

We complete our study by calculating branching ratios and partial decay width ratios 
also for all relevant D — > VlVi decays. They are listed in Table |2] together with known 
experimentally measured values. 

Finally, we summarize our results: We have investigated semileptonic form factors 
for D — > P and D — > V decays within an approach which combines heavy meson and 
chiral symmetry. The form factors are parametrized to satisfy all constraints coming 
from HQET and SCET. The contributions of excited charm meson states are included 
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FIGURE 1. q 2 dependence of the D° -> K (upper left), D° -> %~ (upper right) and D s -»■ K° (lower) 
transition form factors. 



TABLE 2. The branching ratios and partial decay width ratios for the D — > V semileptonic 
decays. Comparison of our model fit with experiment as explained in the text. 



Decay 


3§(Mod.) [%] 


38 (Exp.) [%] 


ryry (Mod.) 


r+/r± (Mod.) 


D ^K* 


2.2 


2.15 ±0.35 [18,] 


1.14 


0.22 


D ^p 


0.20 


0. 194 ± 0.039 ± 0.013 [25J 


1.11 


0.14 


D+ -> K* 


5.6 


5.73 ±0.35 U£] 


1.13 * 


0.22 1 


D + ^po 


0.25 


0.25 ±0.08 OS] 


1.11 


0.14 


D+ ->(» 


0.25 


0.17±0.06±0.01 [25J 


1.10 


0.14 


D s ^<t> 


2.4 


2.0 ±0.5 usj 


1.08 


0.21 




0.22 


1.03 


0.13 



* Exp. 1.13 ±0.08 1 18] 
"f Exp. 0.22 ±0.06 1 18] 



into analysis. The values of form factors at q^ax are calculated using heavy meson chiral 
Lagrangian. The second poles of the F + (q 2 ), V(q 2 ) and Ao(q 2 ) form factors are saturated 
by the presence of the next radial excitations of D*, \ and . The single pole q 2 behavior 

of the A\(q 2 ) form factor is explained by the presence of a single 1 + state relevant to 
each decay, while in A2(q 2 ) in addition to 1 + states one might include their next radial 
excitations. 

The obtained q 2 dependence of the form factors is in good agreement with recent 
experimental results and existing theoretical studies. The calculated branching ratios are 
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FIGURE 2. q 2 dependence of the D° -v K*~ (upper left), D° -> p " (upper right) and D s -> <j) (lower) 
transition form factors. 



close to the experimental ones. We hope that the ongoing experimental studies will help 
to shed more light on the shapes of the D — > P,V form factors. 
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